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ABSTRACT 



' Context. The new VISual and Infrared Telescope for Astronomy (VISTA) has started operations. Over its first five years it will be collecting 
data for six Public Surveys, one of which is the near-infrared YJK^ VISTA survey of the Magellanic Clouds system (VMC). This survey 
comprises the Large Magellanic Cloud (LMC), the Small Magellanic Cloud, the Magellanic Bridge connecting the two galaxies and two fields 
in the Magellanic Stream. 

Aims. This paper provides an overview of the VMC survey strategy and presents first science results. The main goals of the VMC survey 
are the determination of the spatially-resolved star-formation history and the three-dimensional structure of the Magellanic system. The VMC 
survey is therefore designed to reach stars as faint as the oldest main sequence turn-off point and to constrain the mean magnitude of pulsating 
variables such as RR Lyrae stars and Cepheids. This paper focuses on observations of VMC fields in the LMC obtained between November 
2009 and March 2010. These observations correspond to a completeness of 7% of the planned LMC fields. 

Methods. The VMC data are comprised of multi-epoch observations which are executed following specific time constraints. The data were 
reduced using the VISTA Data Flow System pipeline with source catalogues, including astrometric and photometric corrections, produced 
and made available via the VISTA Science Archive. The VMC data will be released to the astronomical community following the European 
Southern Observatory's Public Survey policy. The analysis of the data shows that the sensitivity in each wave band agrees with expectations. 
Uncertainties and completeness of the data are also derived. 

Results. The first science results, aimed at assessing the scientific quality of the VMC data, include an overview of the distribution of stars in 
colour-magnitude and colour-colour diagrams, the detection of planetary nebulae and stellar clusters, and the Ks band light-curves of variable 
stars. 

Conclusions. The VMC survey represents a tremendous improvement, in spatial resolution and sensitivity, on previous panoramic observations 
of the Magellanic system in the near-infrared, providing a powerful complement to deep observations at other wavelengths. 
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1. Introduction 

The cosmological paradigm for the formation and evolution of 
galaxies suggests that large structures formed as a sequence of 
mergers of smaller objects (White & Frenk [T99T]) . The theo- 
retical framework relies on cold dark matter simulations and 
is supported by high redshift observations (York et al. I2000I I 
and by studies of the cosmic microwave background (Spergel 
et al. 12003) . but the major difficulty is to reproduce the bary- 
onic (stars, gas and dust) content of the Universe. Therefore, 
the study of the assembly process of nearby galaxies via re- 
solved stars is a crucial aspect to understand how structures 
in the Universe form and evolve (Tolstoy et al. 120091 ). In par- 
ticular, dwarf irregular galaxies are well suited because their 
low metallicity and high gas content provide information about 
galaxies at an early stage of evolution. The closest prototypes 
of interacting dwarf galaxies that offer an excellent laboratory 
for this near-field cosmology are the Magellanic Clouds (MCs). 

The Magellanic system is located at a distance of ~ 
57 kpc (e.g. Cioni et al. 2000') and comprises: the Large 
Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC), 
the Magellanic Bridge and the Magellanic Stream. The LMC 
is a dwarf irregular galaxy seen nearly face-on (e.g. van der 
Marel & Cioni .2001) and sometimes referred to as a late- 
type spiral galaxy, rich in gas and actively forming stars. The 
SMC is a highly inclined dwarf irregular galaxy also referred 
to as a dwarf spheroidal galaxy (Zaritsky et al. ,2000 ) with less 
active star formation. The LMC is probably just a few kpc 
thick along the line-of-sight, but the SMC has a more complex 
structure that may extend up to 20 kpc along the line-of-sight 
(e.g. Westerlund 119971 Groenewegen 120001 Subramanian & 
Subramaniam |2009al l. There is a bar embedded in each galaxy 
(Subramaniam & Subramanian 2009b. Gonidakis et al. 2009 ). 
The Magellanic system is metal-poor, the metallicity is about 
1/2, 1/4 and 1/10 that of the Sun for the LMC, SMC and the 
Bridge, respectively. The MCs have experienced an extended 
star formation history (e.g. Hill et al. 2000j Zaritsky et al. 120021 
12004 '; Cole et al.'2005'; Pompeia et al. 20081 Gallart et al. l2008l 
Carrera et al. 2008J . 

The dynamical interaction between the MCs may be re- 
sponsible for the various episodes of star formation (Zaritsky 
& Harris ,2004) and for the creation of the Bridge (Irwin et 
al. II985I Gordon et al. 120091 ) which connects the two galax- 
ies and clearly has young stars associated with it (Irwin [1991] 
Battinelli & Demers 1998). The Stream appears (to date) as 
a purely gaseous feature spanning more than 100 deg in the 
Southern sky (Guhathakurta & Reitzel 119981 ). A tidal origin 
of the Stream from the interaction between the LMC and the 
Milky Way (MW) has been ruled out by new proper motion 
measurements (Kallivayalil et al. '2 006al I2006bl ). Alternative 
explanations are: ram pressure (Mastropietro et al. 12005 ) and 
tidal origin from the interaction between the LMC and the SMC 
(Besla et al. lIOTOl) . 

The interaction between the MCs and the MW is repre- 
sentative of the environmental effects that large galaxies with 
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satellites (low-mass dwarf galaxies) experience elsewhere in 
the Universe. This suggests that the MCs may have entered the 
Local Group as part of an association (cf. Tully et al. 120061 
Moss 2006 Knebe et al. 2006) and that, in the future, a minor 
merger between the LMC and the MW may occur (cf. Ibata et 
al. 119941 12003I I. An alternative is that the MCs may be tidal 
dwarfs expelled during a previous merger event involving M3 1 
(Yang & Hammer l20101 l. although this would imply they have 
retained little dark matter from their parent halos (Barnes & 
Hernquist 1992). At present, basic assumptions are being chal- 
lenged: What is the origin of the MCs? Do the MCs consti- 
tute a binary system and, if so, for how long? Have the MCs 
interacted with the MW or are they on their first approach 
(Besla et al. l2007l D'Onghia & Lake [2008b ? How have the star- 
formation histories of the LMC and the SMC been influenced 
by interaction? Does the geometry of the system depend on age 
and metallicity? How do star clusters form and evolve in the 
MCs? What is the fate of the MCs and will they merge with the 
MW? Will the Bridge evolve into a dwarf galaxy (Nishiyama 
et al. 2007 )? Does the LMC have an ordinary bar and how does 
it influence the LMC evolution? Or, is the offset bar a separate 
galaxy being merged into the LMC disc? Does the LMC have 
a metal poor old halo? Why is there a significant difference in 
structure between the gas and stars in the SMC? Does the SMC 
have a bulge? 

To answer all these questions we must resolve the stel- 
lar populations and study them in detail. A fundamental step 
in this direction has been the many panoramic imaging sur- 
veys that have provided multi-wavelength observations of the 
Magellanic system. Except for the dedicated optical MCPS sur- 
vey (Zaritsky et al. 2002, 2004) information about the overall 
population of the MCs has been obtained from surveys with 
different original goals, including microlensing optical surveys 
(e.g. MACHO - Alcock et al. l2000l EROS - Tisserand et al. 
120071 and OGLE - Udalski et al. 1992 ), and infrared sky sur- 
veys (e.g. IRAS - Schwering il989i 2MASS - Skrutskie et 
al. I2056l and DENIS - Cioni et al. l2000l l. A continuation of 
OGLE is still in progress while other large-scale, near-infrared 
(near-IR) surveys (IRSF - Kato et al. l2007i ) offer somewhat 
more sensitive data than 2MASS. The surveys in the mid- 
infrared with the Spitzer Space Telescope (SAGE - Meixner 
et al. i2006. Bolatto et al. 1 2^71 Gordon et al.[2009^, Bonanos et 
al. 1201 Oi l of the central part of the galaxies and an optical sur- 
vey of the outermost regions (Saha et al. 120101) have recently 
been completed. These surveys have provided data covering 
most of the electromagnetic spectrum, but their common depth 
is limited to moderately bright giant stars. The development 
of the VISual and Infrared Survey Telescope for Astronomy 
(VISTA - Emerson & Sutherland 120 101 ) offers a unique oppor- 
tunity to acquire near-IR data of unprecedented sensitivity in 
the Magellanic system. This is the underlying objective of the 
near-IR YJKs VISTA survey of the Magellanic Clouds system 
(VMC). 

This paper is organised as follows. Section|2]introduces the 
VMC survey and then describes the observing strategy and first 
observations. Section [3] describes the data reduction steps for 
producing images and catalogues for individual observations, 
while Sect, m describes the subsequent stages of reduction for 
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deep and linked observations and presents the archival proce- 
dures. Section |5] shows results from the first data. Section |6] 
concludes this study and the Appendix gives the coordinates of 
the VMC fields. 

2. VMC survey 

The VMCQ is a uniform and homogeneous survey of the 
Magellanic system in the near-IR with VISTA. The main pa- 
rameters of the survey are summarised in Tab.[T] It is the result 
of a letter of intent submitted in 2006 and a science and man- 
agement plan approved early in 2008. The main science goals 
of the survey are the determination of the spatially-resolved 
star-formation history (SFH) and the three-dimensional (3D) 
structure of the Magellanic system. VMC observations will 
detect stars encompassing most phases of evolution: main- 
sequence stars, subgiants, upper and lower red giant branch 
(RGB) stars, red clump stars, RR Lyrae and Cepheid variables, 
asymptotic giant branch (AGB) stars, post-AGB stars, plan- 
etary nebulae (PNe), supernova remnants (SNRs), etc. These 
different populations will help assess the evolution of age and 
metallicity within the system. 

The SFH will be recovered from the analysis of colour- 
magnitude diagrams (CMDs) and simulations of the observed 
stellar populations, accounting for foreground stars and ex- 
tinction. Kerber et al. (I2009I I show a preliminary assessment 
of the SFH accuracy that can be expected from VMC data. 
Modelling near-IR colours bears no greater uncertainty than 
optical colours and the near-IR is particularly sensitive to the 
colour of the oldest turn-off" stars. The 3D geometry will be 
derived using different density and distance indicators like 
the luminosity of red clump stars and the Cepheid and RR 
Lyrae period-luminosity relation, and period-luminosity-colour 
and Wesenheit relations. These results will complement those 
based on 2MASS data for the AGB and upper RGB popula- 
tions of the LMC (van der Marer2001) and those from optical 
data of the central regions of the galaxies that are affected by 
a higher reddening (e.g. Subramanian & Subramaniam 120 1 Oi l 
and crowding. We will constrain the epoch of formation of each 
galactic component by mapping the extent of different kinds of 
stars and deriving, using up-to-date stellar evolutionary mod- 
els, ages and metallicities. The VMC survey science addresses 
many other issues in the field of star and galaxy formation 
and evolution, such as: stellar clusters and streams; extended 
sources; proper motions; star formation; distance scale; models 
of Magellanic system evolution; extinction mapping. 

2.1. The yiSTA telescope and camera 

VISTA is a new 4m class telescope developed in the United 
Kingdom (UK). It formed part of the in-kind contribution of 
the UK to joining the European Southern Observatory (ESO). 
The telescope has an alt-azimuth mounting and is located just 
1 .5 km from the Very Large Telescope (VLT) site. The VISTA 
infrared camera (VIRCAM) is equipped with an array of 16 
Raytheon detectors with a mean pixel size of 0.339" and a field 
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of view of 1 .65 deg^. The VIRCAM has a set of broad-band fil- 
ters: Z, Y, J, H and and a narrow-band filter at 1 . 1 8 fim. The 
point spread function of the system is specified to have 50% 
of the light from a point source contained within a circle of di- 
ameter 0.51". The telescope and its camera are described by 
Emerson et al. (2006) and Dalton et al. ('2006'), while the per- 
formance during commissioning is presented by Emerson et al. 
(1201 Oi l. The science verification programmes are summarised 
by Arnaboldi et al. (|20T0] |. 

VISTA is the largest wide-field near-IR imaging telescope 
in the world and it is designed to perform survey observations; 
at present there are six Public Surveys underwajH VISTA ob- 
serves a continuous area of sky by filling in the gaps between 
the detectors using a sequence of six offsets, each by a signif- 
icant fraction of a detector. The combined image corresponds 
to a VISTA tile that covers ~ 1 .5 deg^, while individual offset 
positions are named paw-prints and cover an area of 0.59 deg^. 
The resulting VISTA tile observes each part of sky within it at 
least twice, except for two edge strips in the extreme 'Y' di- 
rections of the array. The tiling process for the VMC survey 
ensures that adjacent tiles overlap sufficiently to provide two 
observations of these areas as well. The combination of offsets 
generates small overlapping areas with exposures as large as 
six times a single paw-print. These steps are described in the 
VISTA user manual 

The VMC survey is a near-IR YJKs survey of the 
Magellanic system. Compared to previous surveys using the 
2MASS and DENIS telescopes, the choice of VISTA filters 
for the VMC survey was driven by the following considera- 
tions. A total of three filters were chosen to be able to analyse 
colour-colour diagrams. The availability of a large colour spac- 
ing to allow for a good characterisation of the subgiant branch 
population to derive the SFH. The relation between J - H and 
7 - A's is quite linear for MC giants, so observations in the 
H band would have provided very little information and prior- 
ity was given to a bluer filter. The J band was preferred over 
H as it also suffers less from atmospheric effects, and provides 
comparative observations with respect to previous surveys. The 
choice of Y over Z was motivated by a reduced confusion limit. 
The Ks filter was required to determine the average magnitude 
of variable stars that, at this wavelength, obey a clear period- 
magnitude relation that is unaffected by other stellar parame- 
ters that degrade the relation at bluer wavelengths (Sect. l572] i. A 
comparison among the filter transmission curves for the VMC, 
2MASS and DENIS observations is shown in Fig. [T] The ex- 
posure time for the VMC survey (Tab. [TJ is designed to meet 
the two key scientific objectives: the SFH and the 3D geometry. 
An accurate determination of the SFH requires CMDs reaching 
the oldest main sequence turn-off, to allow for sampling of dif- 
ferent stellar populations. The investigation of RR Lyrae stars 
and Cepheids requires monitoring observations across specific 
time intervals. The stacking of the observations that are needed 
for deriving the mean magnitude of variable stars in the 
band meets also the depth requirement for the SFH. The split of 
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Filter 


Y 


/ 




Filter 


Y 


J 




Central wavelength O^m) 


1.02 


1.25 


2.15 


Exposure time per epoch (sec) 


800 


800 


750 


Bandwidth (yum) 


0.10 


0.18 


0.30 


Number of epochs 


3 


3 


12 


Detector Integration Time - DIT (sec) 


20 


10 


5 


Total exposure time (sec) 


2400 


2400 


9000 


Number of DITs 


4 


8 


15 


Predicted sensitivity per epoch (Vega mag) 


21.3 


20.8 


18.9 


Number of exposures 


1 


1 


1 


Signal-to-noise per epoch at depth required 


5.7 


5.9 


2.9 


Micro- stepping 


1 


1 


1 


Total predicted sensitivity (Vega mag) 


21.9 


21.4 


20.3 


Number of Jitters 


5 


5 


5 


Total signal-to-noise at depth required 


10 


10 


10 


Paw-prints in tile 


6 


6 


6 


Saturation limit (Vega mag) 


12.9 


12.7 


11.4 


Pixel size (arcsec) 


0.339 


0.339 


0.339 


Area (deg-) 


184 


184 


184 


System FWHM 


0.51 


0.51 


0.51 


Number of tiles 


110 


110 


110 



Jitter pattern = JillerSn. Tile pattern = Tilebzz- 



T-^ I I I r-7Ti I I I I I I I I r 




\{/j.m) 

Fig. 1. Filter transmission curves for the VMC survey {YJK^ - 
black continuous lines) compared with the transmission of the 
2MASS {JHK, - red dashed lines) and DENIS {UK, - blue 
dotted lines) surveys. 



epochs in the Y and J filters is instead purely driven by schedul- 
ing requirements. 

2.2. VMC area coverage 

Observing the entire Magellanic system extending over hun- 
dreds of kpc (the Stream covers half the sky) is a daunting task. 
The VMC survey concentrates, therefore, on a moderately ex- 
tended area (~180 deg^) that includes the classical diameter 
limit at Z? 25 mag arcsec"^ for both galaxies (Bothun & 
Thompson II 988l l as well as major features traced by the dis- 
tribution of stars (e.g. Irwin 119911 Bica et al. 2008i and H I 
gas (e.g. Staveley-Smith et al. 12003! Hatzidimitriou et al. 120051 



Muller et al. I2003I I. the Bridge and two fields at specific loca- 
tions in the Stream. 

The LMC area (116 deg^) is covered by 68 tiles, while 27 
tiles cover the SMC (45 deg^) and 13 cover the Bridge (20 
deg^), see Fig. |2] Additionally, 2 tiles (3 deg^) are positioned 
in the Stream, one approximately to the North of the centre 
of the Bridge (corresponding to a dense area of gas) and the 
other, at similar right ascension, to the North of the SMC, cor- 
responding to a dense area of stars following the simulations by 
Mastropietro (I2009I I. VMC tile centre coordinates are given in 
the Appendix. Each tile is identified by two numbers: the first 
number indicates the row and the second the column of the po- 
sition of the tile in the mosaic that covers the system. Note that 
a separate tiling pattern has been defined for each region (Figs. 
IA.lilA.2llA.3l l. Row numbers increase from South to North and 
column numbers increase from West to East. 

Tiles covering the LMC were oriented at a position angle 
of -1-90 deg. The default orientation (position angle = deg) 
points the 'Y' axis to the North and the 'X' axis to the West. 
The position angle is defined to increase from North to East. 
This represents the best compromise between minimum num- 
ber of tiles and maximum area, increasing the efficiency of the 
survey. The overlap between the doubly-covered sky areas in 
adjacent tiles corresponds to 60" in both 'X' and 'Y' direc- 
tions. The LMC mosaic was created using the Survey Area 
Definition Tool (SADT - Arnaboldi et al. 2008). A geodesic 
rectangle centred at o' = 05:35:50, 5 = -69:27:12 (J2000), 
with width =11.8 deg and height = 15.9 deg, was created as 
the basis of the tiling process. Outer tiles were removed leading 
to the pattern shown in Fig.|2](see also Fig. lA.lb . The area cov- 
ered by the tiles was checked against the distribution of stellar 
associations, carbon stars and other stellar objects using Aladin 
(Bonnarel et al. 1 20001 1. The centre of the rectangle was adjusted 
to include the 30 Doradus nebulosity within a single tile and, 
similarly, the field that the future space mission will re- 

peatedly observe for calibration. 

In the process of creating the mosaic, SADT requires as 
input the observing parameters that are associated to small 
(i.e. jittering) and large (i.e. mosaicking) displacements in the 
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tile position. For the VMC survey the maximum jitter was set 
to 15", the backtrackStep to 100 and the tiling algorithm to 
Tile6zz (these parameters are described in the SADT user man- 
ual). Guide stars were assigned automatically to each tile using 
the GSC-2 reference catalogue (Lasker et al. 12008b . This pro- 
cess may result in shifting the tile centre in case an insufficient 
number of reference stars is available, but this was not the case 
for LMC tiles. 

Tiles covering the SMC region were placed at a position an- 
gle of deg. Keeping the wide tile-edge approximately along 
the right ascension direction produces maximum coverage for 
a minimum number of tiles. This also implies centring the 
geodesic rectangle at a = 00:50:00, 6 = -73:00:00 (J2000), 
with width =8.0 deg and height =8.0 deg. This rectangle rep- 
resents the basis of the tiling process and outer tiles were sub- 
sequently removed leading to the pattern shown in Fig. |2] (see 
also Fig. lA.2l ). As before, the area covered by tiles was checked 
against the distribution of different stellar objects. The posi- 
tion of the centre of the rectangle was tuned to match the area 
that will be observed in the optical domain by the VSllfl (VLT 
Survey Telescope) as part of the STEP survey (P.I. Ripepi; 
Capaccioli et al. 12005b and to provide sufficient overlap for a 
consistent calibration with the Bridge area. 

Previous observations of the Magellanic Bridge by 
Battinelh & Demers (,1992j and Harris (.2007J covered fields 
departing from the LMC and tracing arcs at different declina- 
tions with only a few fields between the two. In order to max- 
imise the total population of stars that the VMC survey will de- 
tect the 2MASS, DENIS and SuperCOSMOsH databases were 
explored for stellar densities and compared with the previous 
observations. VMC tiles were then positioned to overlap with 
the area that provides a good sampling of the stellar population 
while also following a continuous pattern. Following a simi- 
lar procedure for the LMC and SMC areas, a geodesic rectan- 
gle was drawn with centre at a = 03:00:00, 6 = -74:30:00 
(J2000), with width 13.5 deg and height 3 deg. Outer tiles were 
subsequently removed leading to the pattern shown in Fig. [2] 
(see also Fig. lA.3l ). 

The two tiles positioned in the Stream region were prepared 
using the same parameters as for LMC tiles, except for the po- 
sition angle set to deg. Two geodesic rectangles centred on 
the two fields were defined, each with a default size equal to a 
single tile. 

2.3. VMC observations 

The VMC data described here were mostly obtained during 
Science Verification observations (15 - 31 October 2009) and 
the so-called dry-run period (1 November 2009 - 31 March 
2010) when VISTA was tested and survey operations were still 
being defined. A small amount of data was also taken during 
ESO Period 85 (1 April 2010 - 30 September 2010). The bulk 
of the VMC observations will be caiTied out during the odd- 
numbered ESO periods starting in October every year and end- 
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Table 2. VMC required weather conditions. 



Band 


Seeing 


Seeing 




(arcsec) 


(arcsec) 




uncrowded 


crowded 


V 

/ 


1 9 


1 n 

l.U 


J 


1.1 


0.9 




1.0 


0.8 


Field 


Tile row 


Airmass 


LMC 


2,3,4 


1.7 


LMC 


5,6,7 


1.6 


LMC 


8,9,10,11 


1.5 


Stream 


1,2 


1.5 


SMC 


2, 3, 4, 5 


1.7 


SMC 


6,7 


1.6 


Bridge 


1,2,3 


1.7 



Moon distance > 80 deg. Sky transparency = variable, thin cirrus. 

ing in March the following year because of the seasonal ob- 
servability of the Magellanic system. 

Observations of the VISTA Public Survey are obtained in 
service mode by ESO staff. This guarantees efficiency of op- 
erations and a high level of data homogeneity. The requested 
observing conditions for the VMC survey are summarised in 
Table |2l A dozen tiles, centred on the most crowded regions 
i.e. 30 Dor and the central regions of both LMC and SMC, 
have more stringent seeing conditions. This is necessary to pre- 
vent confusion in the bluest bands, while band observations 
will not be limited by confusion for a seeing < 0.9". The best 
FWHMs in the VISTA images are 0.6 - 0.7" and undersam- 
pling, with respect to a pixel size of 0.334", is not a cause of 
concern in the data treatment. 

The Magellanic system never rises above 50° from the hori- 
zon. Therefore, a compromise had to be made between observ- 
ing at reasonable airmass and achieving continuous observabil- 
ity over about five months for the monitoring process. The max- 
imum airmass constraints were optimised as a function of the 
tile declination, as shown in Tab. |2] A violation of any observ- 
ing constraint by 10% is still considered as if observations were 
obtained within specifications. It is also assumed that the air- 
mass constraint may be violated by more than 10% provided 
the seeing constraint is not. 

Table[T]describes the main parameters of the VMC observa- 
tions. The total exposure time is calculated as follows: (number 
of epochs) X 2 X (number of jitters) x (number of DITs) x DIT. 
The factor of two comes from the tiling pattern, during which 
most points of the sky are observed twice (on average). The ex- 
ceptions are the tile edges, observed once, and some areas of 
extra overlap among the detectors that are observed four or six 
times. For example, for the band the total exposure time is: 
12x2x5x 15 x5 = 9000 sec. 

The observations follow the nesting sequence FPJME (see 
VISTA User ManuaJl). This sequence first sets a filter, then 
obtains images at all jittered positions of the first paw-print, 

' http://www.eso.org/sci/facilities/paranal/instruments/vista/doc/ 
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a (J2000) 



Fig. 2. Magellanic system area tiled for VMC observations. The underlying image shows the HI distribution (McClure-Griffiths 
et al. 1200 9). VISTA tiles are colour coded as follows. Blue rectangles represent tiles for which observations have started during 
the dry-runs and P85, green rectangles are for tiles with observations in P86, and red tiles are observations that will not being 
before P87. 



before moving to the next paw-print and taking all the jittered 
images at that position, and so on, until all six paw-prints that 
form the tile are completed. The jittern pattern was jitterSn and 
the tile pattern was Tile6zz- 

The progress of the VMC survey is shown in Fig. |2] and 
detailed in Tab. |3] The reader is referred to the VMC' public 
web page for following up the survey progress beyond that de- 
scribed in this paper. 

2.3.1 . Science Verification 

The observation of a single paw-print in the YJK^ filters was 
requested during the instrument Science Verification (S V) time. 
The main goal was to verify the observing strategy and to test 
the saturation limits for the adopted DITs of 6, 10 and 20 sec 
in the K^, J and Y bands, respectively. The observations were 
carried out on the night of 17 October 2009, but the data were 
only available after the dry-runs had started. 

The SV data indicate that the saturation limit varies from 
one detector to another On average, it is 11 .6 mag for the 6 sec 
Ks band exposure, with seeing of 0.94". The linearity is mild 
at 1 1 .6 mag, but becomes severe for point sources brighter than 
10.7 mag. These limiting magnitudes are somewhat fainter than 
predicted. Prior to the beginning of the dry-runs it was decided 
to reduce the band DIT down to 5 sec, which improves the 
photometry of the bright stars and maximises the range of over- 
lapping magnitudes with 2MASS. The observing strategy was 



Table 3. VMC survey progress. 



Description 


VMC 


LMC 


LMC 




all 


all 


season I 


Total number of tiles 


110 


68 


6 


Total number of epochs 


1980 


1224 


108 


Total number of F-band epochs 


330 


204 


18 


Total number of y-band epochs 


330 


204 


18 


Total number of A^s-band epochs 


1320 


816 


72 


Observed number of F-band epochs 


18 


18 


18 


Observed number of y-band epochs 


18 


18 


18 


Observed number of A^s-band epochs 


51 


51 


51 


Total number of observed epochs 


87 


87 


87 


Completion in the F-band 


5.4% 


8.8% 


100% 


Completion in the y-band 


5.4% 


8.8% 


100% 


Completion in the /fs-band 


3.9% 


6.2% 


71% 


Total completion 


4.4% 


7.1% 


80% 



then established and no modification was needed for the obser- 
vations that had akeady started. 

Other SV observations were obtained on 28 Nov 2009 in 
the SMC. These are band observations of one tile, with a six 
paw-print mosaic and DIT= 10 sec, with the purpose of check- 
ing the sky-subtraction procedures. Images were obtained to 
test an off-sky algorithm which is not applicable to VMC data. 
In addition, some 2MASS touchstone fields that are observed 
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for photometric checks are located in the VMC survey area, 
providing extra, albeit shallow (DIT= 5 sec), multi-epoch data 
in all VIRCAM filters. 

2.3.2. Dry-runs during 1 November 2009 - 14 
February 2010 

One hundred hours of VMC observations were submitted for 
execution during the 1 November 2009-14February2010dry- 
run period and 70% were executed. Tables |4] and |5] summarise 
the observations where each date corresponds to the observ- 
ing date of one observing block (OB). Each VMC OB includes 
only one filter. These observations were prepared using a new 
version of the Phase II Proposal Preparation^ tool especially re- 
vised for Public Surveys with the VISTA and VST telescopes. 
It allows the user to have more control over the survey exe- 
cution with new high-level tools called scheduling containers 
(Arnaboldi et al. '2008). Three types of containers were avail- 
able: concatenations (grouping together OBs for back-to-back 
execution), time-links (imposing time constraints for the exe- 
cution of OBs, including execution of OBs in a user-defined 
sequence), and groups (that improve the prioritisation of OBs 
to ensure the completion of one set of OBs before another set 
is started). 

Six VMC fields were observed in the LMC (Fig. |2]i. One 
field covers the famous 30 Dor region (6_6), one field corre- 
sponds to the South Ecliptic Pole (SEP) region that Gaia will 
observe repeatedly for calibration (8_8), a pair of fields are lo- 
cated in the Northern outer part of the LMC disc (8_3 and 93, 
overlapping in declination) while the remaining two fields over- 
lap in right ascension and are located towards the Bridge (4_2 
and 4_3). The 30 Dor field was chosen because it is well stud- 
ied and represents a crowded region of the LMC rich in stars 
and gas. The Gaia field was chosen for its importance to the 
Gaia mission; early analysis and release of these data will ben- 
efit the astronomical community in general. The two fields in 
the Northern disc were chosen to sample an uncrowded and 
external area of the galaxy. Those towards the Bridge satisfy 
a similar criterion but sample a different region of the LMC 
disc. Overall the six fields were also chosen to reflect different 
coordinates and observing conditions. 

Table |4] shows the observations of concatenations and 
groups for each LMC field. Groups represent observations in 
the Y and J bands that are not linked by time, i.e. for a given 
field they can be executed on any night including on the same 
night. They are labelled Yl, 71, Y2 and 72. Concatenated ob- 
servations are preceded by the letter C and the integration time 
per filter is shorter than for a single VMC epoch at that fil- 
ter (Tab. [TJ, but the different filters are observed back-to-back. 
The time Umit for a concatenation, during this period, was 1 .5 
hours. Two concatenations were, therefore, necessary for re- 
covering one VMC epoch per filter. An important feature of 
concatenations is that the sky condition requirements are met 
during its entire duration. This is not always the case and then 
the concatenation has to be repeated. This explains why in Tab. 
|4]there are additional observations of OBs that are part of con- 



catenations but without observations in all three filters. These 
additional observations will be used to investigate source vari- 
ability and to produce deep stacks (Sect. [Hi if the observing 
condition criteria (Tab. |2| were met. 

Table |5] shows time-link sequences where 1 1 individual 
band OBs are associated with each LMC field. Each sequence 
can start at any time, independent of the execution of observa- 
tions listed in Table [H i.e. a TKl OB can be observed on the 
same night as an OB in the Y, J or A's filter taken as part of a 
concatenation or group. Once started, the next observation of a 
TKn OB, in a time-link sequence, is obtained at intervals equal 
or larger than: 1, 3, 5 and 7 days for epochs 2 to 5, respectively, 
and thereafter at least 17 days from each previous observation 
(epochs 6 to 11). Some OBs in the time-link sequence were 
also repeated. 

2.3.3. Dry-runs during 15 February 2010-31 IVIarcli 
2010 and ESO Period 85 

The remaining dry-run period until the end of March 2010 
and the subsequent ESO survey period (P85) until the end of 
September 2010 were not ideally suited for observations of 
the Magellanic system, except for short windows at the begin- 
ning and the end of the period. The VMC survey was assigned 
30 hours during this period. Seven hours were scheduled to 
progress on the Ks monitoring of four of the six LMC fields 
started during the first dry-run period. The remaining hours 
were for observations of the fields in the Stream. The over- 
all observation of the LMC fields during the dry-runs and P85 
bring the completion rate of the six fields to 80%. 

The strategy for the fields in the Stream was changed be- 
cause concatenations were limited to a total execution of time 
of 1 hour. This limit imposes a further subdivision of the VMC 
epochs that may increase the observing overheads. We have 
decided to adopt a different procedure that keeps the overheads 
more or less constant and does not interfere with the require- 
ment of observing different filters as close in time as possible. 
In the new approach there are 3 concatenations per field and 
each one contains only two filters: YJ, YK^ and JK^. The sum 
of the integration times per filter is the same as that from the 
two previous concatenations and it is also the same as a VMC 
epoch. The first observations of the Stream will be described in 
a subsequent paper of this series. 

3. Data Reduction 

The raw VISTA images acquired for the VMC survey were 
reduced by the VISTA Data Flow System (VDFS) pipeline 
at the Cambridge Astronomical Survey Unit (CASLEI). The 
VDFS pipeline is specifically designed for reducing VISTA 
data (Irwin et al. 120041 ) and is used to process up to 250 
GB/night of data. The pipeline is a modular design allow- 
ing straightforward addition or removal of processing stages. 
The VMC data are reduced together with other VISTA data 
on a weekly basis. Prior to this science reduction the data are 
checked at the observatory site (ESO Chile) using a simplified 



http://www.eso.org/sci/observing/phase2/P2PP/P2PPSurveys.html 
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Table 4. LMC epochs; concatenations (C) and YJ bands. 



Tile 


CKl 


CJ\ 


CYl 


CKl 


CJ2 


CYl 


Y\ 


71 


Y2 


J2 


4.2 


20.12.09 


20.12.09* 


20.12.09* 


15.01.10 


15.01.10 


15.01.10 


15.12.09 


14.11.09 


16.12.09 


28.11.09 










4.01.10 


4.01.10 


4.01.10 






15.12.09 




4.3 


23.12.09 


23.12.09" 


23.12.09 


16.01.10* 


16.01.10 


16.01.10 


9.12.09 


18.11.09 


10.12.09 


28.11.09 




22.12.09 


22.12.09 


22.12.09 
















6.6 


5.11.09 


5.11.09 


5.11.09 


20.11.09 


20.11.09 


20.11.09 


4.11.09 


4.11.09 


8.11.09 


8.11.09 










26.11.09 


9.11.09 












8.3 


26.11.09 


26.11.09 


26.11.09 


22.11.09 


22.11.09 


22.11.09 


19.11.09 


15.11.09 


21.11.09 


26.11.09 


8.8 


21.11.09 


21.11.09 


21.11.09 


20.11.09 


20.11.09 


20.11.09 


5.11.09 


9.11.09 


9.11.09 


11.11.09 




26.11.09 


26.11.09 












6.11.09* 


7.11.09* 




9.3 


21.11.09 


21.11.09 


21.11.09 


22.11.09 


22.11.09 


22.11.09* 


30.11.09 


1.12.09 


1.12.09 


2.12.09 



" reduced number of jitters and/or paw-prints 
* seeing and/or ellipticity too high 

Table 5. LMC epochs: time linked (T) A's-band monitoring. 



Tile 



TKl TK2 TK3 TK4 TK5 TK6 TKl TKS TK9 TKW TKll 



4.2 
4.3 
6.6 

8.3 



9.3 



14.12.09 17.12.09 5.01.10 14.01.10 22.01.10 
10.12.09 17.12.09 21.12.09 27.12.09 18.01.10 

8.11.09 12.11.09 17.11.09 29.11.09 7.12.09 26.12.09 13.01.10" 31.01.10 19.02.10 11.03.10 10.11.10 
6.11.09" 

3.12.09 6.12.09 22.12.09 28.12.09 14.01.10 31.01.10 22.02.10 14.11.10* 

14.11.09 19.11.09 25.11.09 30.11.09 7.12.09* 25.12.09 14.01.10 31.01.10 19.02.10 08.03.10 26.11.10 

16.11.10* 

4.12.09 9.12.09 19.12.09 23.01.10 31.01.10 24.02.10 

16.01.10* 20.02.10* 



" reduced number of jitters and/or paw-prints 
* seeing and/or ellipticity too high 



version of the VDFS pipeline and library calibrations. The data 
are subsequently checked at ESO in Garching for monitoring 
the instrument performance and to feed updated information 
back to the observatory. 

The most relevant VDFS steps for the reduction of VMC 
survey data are as follows: 

- reset, dark, linearity, flat and de-striping (removing hori- 
zontal stripes in the background) correction; 

- sky background correction (tracking and homogenisation 
during image stacking and mosaicking); 

- jittered and paw-print stacking; 

- point source extraction; 

- astrometric and photometric calibration (the latter put in an 
internally uniform system); 

- bad pixel handling, propagation of uncertainties and eflTec- 
tive exposure times by use of confidence maps; 

- nightly extinction measurements. 

The different observational uncertainties are propagated 
during the data processing, to give the users a clear picture of 
the final data quality. Various quality control parameters are 
calculated during the data reduction to monitor the data, and 
to evaluate both the observing conditions (in retrospect) during 
the observation, and the individual data reduction steps. Among 



them are: the zero-point to measure the atmospheric extinction, 
the FWHM to measure the seeing, the ellipticity to evaluate the 
quality of the guiding and active optics correction, the sky level 
to estimate the background level and its variations, etc.. The 
processing history is recorded directly in FITS headers. 

A tile image is produced by combining 96 different images 
(16 detector images per each of 6 paw-prints). Their sky level 
and individual paw-print astrometric and photometric distor- 
tion are adjusted in the drizzling (combination) process. Tile 
catalogues are produced following the application of a nebu- 
losity filter to the paw-prints in order to remove diffuse vary- 
ing background on scales of 30" or larger (Irwin l20T0b . This 
method has shown that the detection of objects and their char- 
acterisation (astrometry, photometry and morphological classi- 
fication) are considerably improved. 

The first VMC reduced images and catalogues, correspond- 
ing to observations obtained during November 2009, were re- 
ceived in January 2010. These data were reduced with version 
0.6 (vO.6) of the VDFS pipeline and refer to individual paw- 
prints, not yet combined into tiles. This means that in practise, 
for example, the observation of a F-band OB of a given field 
has 6 associated images, 6 catalogues and 6 confidence maps. 
Each image and catalogue are delivered in Rice compressed 
format and are multi-extension FITS files containing each the 
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information for all of the 16 detectors covering the VIRCAM 
field-of-view. In March 2010 we received more data processed 
with an upgraded version of the pipeline (vO.8). These data in- 
clude all observations obtained until the end of January 2010. 
The processed observations obtained after this date and until 
the end of the Magellanic season (March 2010) were received 
in April 2010. Following quality inspection a few images were 
re-reduced (vO.9) to fix some specific problems with a small 
subset of the data. The vO.8 data were ingested into the VISTA 
Science Archive (VSA) by June 2010 (Sect. Ell. They include 
astrometry and photometry in single-band, band-merged and 
epoch-merged tables as well as deep stacks. In October 2010 
vl.O data have become available and the main difference from 
previous releases is that they include tile images and cata- 
logues. 

Figure[3]shows most of a tile of the 6_6 LMC field including 
the star forming region 30 Dor This image was produced for 
an ESO press releas43 where other zoomed-in images are also 
available. The exposure time in the three wave bands was 2400 
sec in Y, 2800 sec in J and 4850 sec in K^. 

3.1. Linearity and background 

The individual VISTA detectors have different non-linearity 
and saturation properties and these properties may also vary 
across a given detector The detector system is non-linear and 
linearity corrections are applied at the pixel level during the 
initial image processing stages. The saturation levels are stored 
in the image header keywords and, together with the peak flux 
derived from the photometry of the observed sources, it is pos- 
sible to establish a correction that enables recovery of stars up 
to a few magnitudes brighter than the saturation limit (Irwin 
I2009I I. This correction is not implemented by VDFS prior to 
the production of VMC catalogues but it is applied at a further 
stage of the data processing (Sect. |4|i. The expected saturation 
values listed in Tab.[T]are sensitive to seeing variations. In fact, 
the 30 Dor region was observed under very good conditions, 
and saturation appears at fainter magnitudes. 

The application of a nebulosity filter to the paw-print im- 
ages prior to the construction of a tile image (Irwin 2010), may 
influence the recovery of the magnitude of stars close to the 
saturation limit. This would, however, only affect heavily satu- 
rated stars, i.e. those where a significant halo with a diameter 
comparable to that of the filter size (~ 30") is visible. 

The sky background for VMC observations is estimated, 
for each paw-print, from all paw-prints observed for one tile in 
a given band and at a given time. This method, referred to as 
'tilesky', has shown very good results even for the 30 Dor tile 
where there is a substantial emission from the nebula. 

Persistence effects due to bright stars are usually automat- 
ically removed by the VDFS pipeline when the observations, 
Hke for VMC, follow the FPJME sequence (see VISTA User 
Manual^). Adding up all VMC images for a given field and 
filter does not produce any noticeable effect due to persistence. 

The moon has a negligible effect on the VMC background 
because it is always >80 deg away from any of the fields. This 

http://www.eso.org/public/news/esol033/ 



results in low contamination, even in the Y band which is most 
susceptible to the lunar contamination. The major absorption is 
caused by water vapour and carbon dioxide in the atmosphere. 
At the VMC filters the background will also be dominated by 
non-thermal aurora emission, OH and O2 lines especially in 
observations obtained 1 .5 - 2 hours after twilight. 

3.2. Astrometry 

Astrometry is based on positions of the many 2MASS sources 
within each detector. The astrometric calibration of a paw- 
print is encoded in the FITS image headers using the Zenithal 
Polynomial projection (ZPN) while a tile refers to a single tan- 
gent plane World Coordinate System (WCS) image (Calabretta 
& Greisen 2002). The median astrometric root-mean-square is 
80 mas and is dominated by 2MASS uncertainties. Residual 
systematic distortions across the VISTA field-of-view are at 
the ~25 mas level and can be further improved, if required, 
by directly characterising the distortion pattern. In a dithered 
sequence the detectors are rotated slightly to maintain the posi- 
tion angle on the sky; a comparison between identical VISTA 
observations shows a residual rotation of ~0.5 pix. 

Figure |4] shows a comparison between the right ascension 
(a) and declination (6) coordinates of stars that are in common 
between VISTA and 2MASS in the VMC field 8.8. The excel- 
lent match shows the quality of the astrometry. A systematic 
shift is perhaps present at the level of ~ 0.01" in both axes. 
This accuracy is perfectly adequate for cross-correlation stud- 
ies between external catalogues and VMC. The relative accu- 
racy within VMC data is higher and a more detailed investiga- 
tion of the suitability of VMC data for studies of, for example, 
proper motions will be addressed elsewhere. 

3.3. Pliotometry 

The photometric calibration relies on the observation of stars 
from the 2MASS catalogue with magnitudes in the range 12 - 
14 mag in all bands. The procedure is similar to that adopted 
to calibrate data from the Wide-Field Camera (WFCAM) 
mounted at the United Kingdom Infrared Telescope (UKIRT). 
For the WFCAM filters, which are very similar to VISTA's (ex- 
cept that WFCAM has K and VISTA has K,), Hodgkin et al. 
(2009) have shown that the calibration of the Y band, not in- 
cluded in 2MASS, is possible where the extinction is not too 
high, i.e. E{B - V) < 1.5. This is well within the average ex- 
tinction values towards the Magellanic system (Westerlund et 
al. I1997I I. However, in star-forming regions the extinction can 
reach larger values and the calibration may not be reliable. To 
remedy this situation a calibration based on the observations of 
standard stars will be produced together with a thorough inves- 
tigation of star-forming regions from previous data. 

A high quality global photometric calibration of the VMC 
survey will be supported by the homogeneity and accuracy of 
the 2MASS catalogue. The best absolute photometry is ex- 
pected to be accurate to about 1% and, on average, 2%, but 
relative photometry will reach a much greater accuracy (~milli- 
magnitudes for brighter sources). At periodic intervals and at 
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Fig. 3. Most of VMC tile of LMC field 6_6. This is a colour composite image where Y is shown in blue, J in green and in 
red. East is to the left and North at the top. The 30 Doradus star forming region is visible together with other smaller regions 
towards the South as well as stellar clusters and the field population. For more details and a high resolution image refer to 
|http://www.eso.org/pubUc/news/esol033/; 
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Fig. 4. Comparison between VISTA and 2MASS astrometry in 
the 8_8 field. Histograms have bins of 0.01" in size. The best 
fitting gaussians are indicated and coiTespond to cr of 0.080" 
and 0.085" for Aa and A5, respectively. 

the end of the survey the global photometric calibration will be 
assessed. 

Figure |5] shows the behaviour of VISTA photometric un- 
certainties in the VMC 8_8 field where data represent stacked 
paw -prints and tiles. Note that the uncertainties reduce by about 
~ 50% compared to the individual tiles, and will reduce fur- 
ther for deep tiles. In the VMC catalogues several aperture flux 
magnitudes are given that sample the curve of growth of all 
images. The recommended aperture 'aper3', used in this work, 
coiTesponds to a core radius of 1" (3 pixels) that contains 75% 
of the total stellar flux in a 0.8" seeing observation. 

Figure |6] shows a comparison between VMC and 2MASS' 
magnitudes. By selecting a suitable range of magnitudes the 
parameters of the comparison (mean, median and sigma) are 
indicated in Tab. |6l Note that by 2MASS* we do not refer to 
2MASS magnitudes, but to the magnitudes obtained using the 
following colour equations: 

J'aMASS- = Jim + 0.550 X (7211, - Ii2m) (1) 
^2MASS- = Jim - 0.070 X (72m - ^2^) (2) 
^2MASS. = + 0.020 X {J2m ' ^2m) (3) 

where 72ni, Him ^2m are 2MASS magnitudes. These are the for- 
mulas used to calibrate the VISTA photometr}{3 They include 
the correct colour term but do not include the small zero-point 
shifts that are necessary to bring the observations into a Vega 
magnitude system. 

" http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/vista- 
sensitivity 



Fig. 5. Photometric uncertainties in the VMC data for stacked 
paw-prints (dashed, dotted and continuous lines in the K^, J 
and Y bands respectively) and tiles (red, green and blue cir- 
cles in the K^, J and Y bands respectively) in the 8_8 field. 
Uncertainties are progressively smaller from to Y and are 
systematically smaller in tiles than in stacked paw-prints. 

Table 6. VMC-2MASS* comparison. 



Band 


Range 


Mean 


Median 


Sigma 




(mag) 


(mag) 


(mag) 


(mag) 


Y 


10.5 - 15.0 


-1-0.018 


+0.012 


0.110 


J 


10.5 - 14.5 


-0.005 


-0.003 


0.092 




10.0 - 14.0 


-0.003 


-0.003 


0.069 



The morphological classification is encoded in a flag with 
values: -1 for stellar, +1 for non-stellar, for noise and -2 
for borderline stellar sources. This classification indicates the 
most probable source morphology and is affected by crowding 
and sensitivity. For details about the photometric apertures, the 
source classification and other catalogue parameters the reader 
is referred to the CASU web pages. 

3.4. Image quality 

The quality of the VMC images is evaluated at three differ- 
ent steps - during the observations, the data reduction, and the 
archiving process. The VMC team performs additional quality 
control checks that interleave with each of these three steps. 
The first identifies obvious causes for re-observation, such 
as observations that exceeded the required constraints or that 
failed to be completed because of technical reasons. These ob- 
servations are usually repeated almost immediately (Tabs.|4l|5]l. 
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2MASS' 

Fig. 6. Magnitude difference between VMC and 2MASS* 
sources for stars in the 8_8 field. Horizontal lines were defined 
within a specific range of magnitudes that is given in the text 
and is indicated in red in the colour version of the figure. 

All images are processed by CASU and archived at the VSA, 
regardless of whether their observing constraints are met. 

All CASU-reduced VMC survey images for individual de- 
tectors have been inspected visually for quality control. The 
purpose of this inspection is to recognise artefacts and residu- 
als from the reduction process but also to identify features that 
are intrinsic to the observations. The results of the quality in- 
spection are as follows: 

- The upper 1/3 of detector #16 is effectively noisy in the 
bluer bands, this also causes a calibration problem for that 
region: so the VMC Y band observations suffer more than 
the J and Kf^ ones. This problem causes an increasing back- 
ground level and influences the detection capabilities. 

- The observations obtained before 20 November 2009 suf- 
fer from intermittent problems in detector #6 channel #14 
(each detector has 16 channels, processed by different 
analog-to-digital convertors) that required replacement of 
a video board in the controller. These stripes cover an area 
comparable to that of a bright foreground star. 

- A special sky frame needs to be used for reducing images 
obtained during the night of 19 November 2009, due to in- 
vestigations on the component responsible for the previous 
problem. This step is implemented in the pipeline reduction 
from version 0.9 onwards. 

- A low quality region at the bottom ('-Y') of detector #4 
creates a horizontal pattern that does not cancel out with 
stacking images obtained from the exposure and jittering 
sequence. This problem will likely not affect subsequent 
reductions of the data. 

- Overall the reduced images show a smooth gradient most 
noticeable in the band (possibly caused by the baffling 



Table 7. Average VMC parameters from all single tile images. 



Band 


FWHM 


EUipticity 


Zero - point 


5(T Mag. Limit 




(arcsec) 




(mag) 


(mag) 


Y 


1.03 ±0.13 


0.065 ±0.011 


23.520 ±0.070 


21.111 ±0.395 


J 


1.00 ±0.10 


0.064 ±0.011 


23.702 ± 0.206 


20.527 ± 0.382 




0.93 ± 0.08 


0.051 ±0.009 


22.978 ± 0.245 


19.220 ±0.340 



system of VIRCAM as a result of thermal radiation), but 
does not present a problem for the source extraction. 

None of these problems require re-observation of the tiles that 

have been obtained for the VMC survey. 

The quality of the images was further inspected by com- 
paring the FWHM of the extracted sources with the expected 
seeing requirements. As expected the FWHM varies among the 
VIRCAM detectors. For example, the average FWHM in the 
corners of detector #1 is larger than in other detectors. This 
implies that some detectors will have a FWHM that exceeds 
the seeing requirement by more than 10%. The average seeing 
among the detectors and the average seeing among the six paw- 
prints of a tile is, however, always within the required Umit, 
except in a few cases, see Tabs. |4] and |5] 

Tables |4] and |5] also show a few cases where an insufficient 
number of paw-prints or jitters are present and tiles are not 
fully sampled. By discarding these problematic observations 
among those that remain only very few need to be re-obtained. 
Table Q shows the average seeing, ellipticity, zero-point and 
limiting magnitude from all VMC images. These parameters 
have been calculated excluding problematic observations and 
represent average values regardless of their execution as part 
of a concatenation, group or monitoring sequence, as well as 
from the specific requirements of the different VMC fields 
with respect to crowding. 



4. Data archive 

The data reduced by the VDFS pipeline at CASU are ingested 
into the VSaEI at the Wide Field Astronomy Unit (WFAU) in 
Edinburgh which is similar to the WFCAM Science Archive 
(Hambly et al. l2008l l. At present, these are data reduced with 
V 1.0 of the CASU pipeline and include all VMC data observed 
until end of May 2010 (Tabs. Hand [S]). At VSA the data are 
curated to produce standardised data products. The software 
that runs at WFAU and populates the VSA is the same that 
runs at CASU and this guarantees that the data are processed 
homogeneously throughout the entire processing chain. 

The most important processes, available at present, for the 
VMC survey are: individual passband frame association and 
source association to provide multi-colour, multi-epoch source 
lists; cross-association with external catalogues (list-driven 
matched photometry); deeper stacking in specified fields; qual- 
ity control procedures. 

http://horus.roe.ac.uk/vsa/login.html 
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There are three main types of VSA tables that are impor- 
tant for the VMC survey. These are the vmcDetection table, the 
vmcSource table and the vmcSynoptic table(s). The vmcDetec- 
tion table contains the catalogues corresponding to individual 
observations. At the moment there is one catalogue per paw- 
print, regardless of band and tile of origin. The vmcSource ta- 
ble contains the list of sources obtained from deep stack im- 
ages and each source is matched in the three VMC bands. In 
practise, each row of the vmcSource table will contain Y, J 
and magnitudes for a source. At present, because VSA is 
organised by paw-prints, the same source may appear two or 
more times in the vmcSource table depending on its location 
with respect to the overlap among the six paw-prints forming a 
tile. The synoptic tables contain the colour information and the 
multi-epoch information for individual observations (a single 
OB). More details about the synoptic tables are given by Cross 
et al. ( I2009I I. The position and magnitude for each source in 
any given table refers to the astrometrically and photometri- 
cally calibrated measurements using the parameters specified 
in the image headers. These parameters are discussed in Sect. 
[3] In addition, there are several quality flags that are specifically 
introduced at the VSA level. These flags identify problems oc- 
curring during the ingestion of the data into the archive, incom- 
pleteness in the set of data (for example missing exposures in 
a paw-print sequence), problems related to the pairing of data, 
etc. 

The magnitudes of the brightest stars are corrected for sat- 
uration effects (Sect. 13. 11 1. Figure [T] shows the magnitude 
difference for stars in the 8_8 field compared to 2MASS before 
and after saturation correction. Saturation effects are present for 
A's < 12 mag, with the magnitudes of brighter stars recovered 
to at least - \Q mag after correction. 

The VSA is queried using Structured Query Language 
(SQL) and a point-and-click web form for browsing. This is a 
dual (sophisticated and simple) end-user interface for the data. 
A key feature to note is the design with multi-waveband cata- 
logue data that allows the user to track back to the individual 
source images and merged-source tables, and present the user 
with a generally applicable, science-ready dataset. The VSA 
has a high-speed query interface, links to analysis tools such as 
TopCat, and advanced new VO services such as MySpace. The 
VSA supports a different range of queries and the most com- 
mon for the VMC survey are: (i) querying the archive to check 
which data have been ingested, (ii) querying the vmcSource ta- 
ble to extract magnitudes from the deep stacks, (iii) querying 
the synoptic tables to extract light-curves and statistics on the 
levels of variability and (iv) querying the VSA using an input 
list of sources and searching for their VMC counterpart. For 
each source in (i) it is possible to inspect postage stamp images 
in each wave band. 

Before creating the deep stack images, the quality control 
results discussed in Sect. l3.4l need to be taken into account. In 
particular, all images except those with a large seeing and/or el- 
lipticity listed in Tabs. |4] and |5] are included in the deep stacks. 
This means that tile images with a reduced number of jitters or 
paw-prints are included into the deep stacks if their observing 
conditions are met. The latter will still be available as individ- 
ual epochs, and will be linked to the other observations in the 




10 12 14 

K 

Fig. 7. Magnitude diff'erence for stars in the 8_8 field in com- 
mon between VMC and 2MASS before (bottom) and after 
(middle) correcting the magnitude of VMC stars approaching 
the saturation limit, and (top) the correction itself developed by 
Irwin (2009). 

synoptic tables, because they may contain useful information 
for variable stars or for source confirmation. 

In the future it may be possible to automatically join the 
under-exposed areas at the two ends of each tile with those of 
the adjacent tiles prior to source extraction. The VSA also con- 
tains external catalogues, like 2MASS, that can be linked with 
the VMC data via an SQL query. Catalogues that are specifi- 
cally important for the Magellanic system, e.g. the MCPS and 
the SAGE catalogues are also being ingested into the VSA. 
VMC is intrinsically a multi-wavelength project and a large 
fraction of its science will come from the linking of VISTA 
data with those from other surveys; the VSA is designed to en- 
able such links. 

5. Analysis and Results 

FigurelHshows the CMDs of the VMC data in the 6.6 (30 Dor) 
and 8_8 (Gaia) LMC fields. These data were extracted from the 
VSA. The magnitudes and colours of each source correspond to 
a single detector, i.e. if the same source was detected in another 
detector it is not included in the CMDs. This is because at this 
stage the archiving process is organised by paw-prints. When 
tiles, resulting from the combination of six paw-prints, become 
available at VSA then the exposure time per source will be at 
least doubled. At present the exposure times per band for the 
sources shown in the CMDs correspond to 1200 sec in Y, 1400 
sec in J and to ~ 4000 sec in Ks- 

The distribution of stars in the CMDs shows clearly the dif- 
ferent stellar populations characterising these LMC fields. The 
blue-most conic structure bending to red colours at bright mag- 
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Fig. 8. Colour-magnitude diagrams of VMC sources in part of the 6_6 (30 Dor) and 8_8 (Gaia SEP) LMC fields. 



nitudes is formed by main-sequence (MS) stars of increasing 
mass with increasing brightness. The MS joins, via the sub- 
giant branch, the RGB beginning at ~2 mag below the red 
clump, the approximately circular region described by the high- 
est concentration of stars. The structure of the red clump de- 
pends on stellar parameters (age and metallicity) but also on 
interstellar extinction. Extinction causes the clump to elongate 
to red colours, as seen in the CMDs for the 30 Dor field where 
extinction is higher than in the Gaia field. The RGB contin- 
ues beyond the red clump at brighter magnitudes describing 
a naiTow structure bending to red colours. The abrupt change 
in source density at the tip of the RGB marks the transition 



to brighter AGB stars. The broad vertical distribution of stars 
below the RGB is populated by MW stars. In the CMDs of 
the Gaia field these are easily distinguished from LMC stars. 
Cepheid and supergiant stars occupy the region of the diagram 
to the bright and blue side of the RGB while RR Lyrae stars 
are somewhat fainter than the red clump and lie more or less 
parallel to the sub-giant branch. 

Figure|9]shows the colour-colour diagram of the VMC data 
in the 8_8 field. The data shown here are the same as in the 
CMDs (Fig. [8]l described above. The distribution of sources 
in the colour-colour diagram resembles the body of an ant. 
Following this analogy: 
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CMDJ-Ks vs. Ks 




LJ '. I I L 

Fig. 9. Ant diagram: the colour-colour diagram of VMC 12 

sources in part of the 8_8 LMC field. (J-Ks) 



- The lower part of the ant body (gaster) at -0. 1 < (Y - J) < 
0.2, -0.3 < (J - Ks) < 0.2, corresponds to the location of 
MS stars in the LMC, with the youngest stars being at the 
bluest extremity. The width of this feature is mainly caused 
by photometric uncertainties. Its extension to the red is also 
set by the limited depth of the VMC data, since the faint MS 
should continue to even redder colours. 

- The middle part of the ant body imesosoma) at 0.2 < {Y - 
J) < 0.4, 0.3 < (J - Ks) < 0.65, coiTesponds to the main 
locus of helium-burning giants in the LMC. The bulk of 
them are in the red clump (see Fig. [8]), but also brighter 
helium burning giants, and stars in the faint extension of 
the red clump, fall in this same blob. 

- The petiole is a small thin extension of the mesosoma at its 
red side, at 0.2 < (Y - J) < 0.3, 0.2 < (J - Ks) < 0.3, and 
is mainly caused by bright stars in the MW foreground - 
more specifically by the intermediate-age and old turn-offs 
of MW disk populations (at Ks < 15). 

In addition to this well-defined petiole, gaster and meso- 
soma are also connected by the relatively less populated 
LMC sub-giants, lowest-luminosity RGB stars, and horizontal 
branch stars in the LMC. 

- The upper part of the ant body (head) is a more complex 
feature. Its main blob at 0.4 < (Y-J) < 0.6, 0.6 < ( J-Ks) < 
0.8 is defined by low-mass stars in the MW foreground, 
especially those with masses < 0.5 Mq which clump at 
the same near-IR colours (J-Ks ~ 0.7; see Nikolaev & 
Weinberg ,2000. Marigo et al. I2003I I. The same structure 
forms a marked vertical feature in the CMDs. 

- Two antennae depart from this head, the upper one at 0.5 < 
(Y-J) < 0.6, 0.8 < (J - Ks) < I being formed by the 
more luminous RGB stars in the LMC, close to their tip of 



Fig. 10. CMD for a region of 4000 x 4000 pixels extracted from 
the 8_8 LMC deep tile. Objects with J-Ks> I and fainter than 
Ks - 16 mag are background galaxies. 

the RGB, extending up to (J - Ks) - 1 mag. This upper 
antenna finishes abruptly because the tip has been reached. 
The lower antenna at (J - Ks) ~ 0.8, (Y-J) > 0.6, is 
more fuzzy, and corresponds to the Y-J red-ward extension 
of low-mass stars in the MW foreground. This red ward 
extension is partially caused by photometric uncertainties 
and by the particular colour-colour relation followed by the 
coolest M dwarfs. 

5.1. Completeness 

Some key science goals of the VMC survey require accurate 
estimates of the completeness of the stellar photometry as a 
function of location across the Magellanic system, and posi- 
tion in the CMDs. This is estimated via the usual procedure 
of adding artificial stars of known magnitudes and positions to 
the images, then looking for them in the derived photometric 
catalogues. 

For this work, the paw-print images, available for the 8_8 
LMC field, were combined to produce a tile image using the 
SWARP tool (Benin et al. I2002I I. A region of the tile image 
with a size of 4000 x 4000 pixels was selected. Then, PSF 
photometry using the DAOPHOT and ADDSTARS packages 
in IRAF was performed. Figure [10] shows an example of the 
photometry on the (preliminary) stacked image. 

The artificial stars are positioned at the same random po- 
sition on the YJKs images. Their mutual distances are never 
smaller than 30 pixels, so that the process of adding stars does 
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Fig. 11. Completeness results for the 8_8 LMC tile with a single 
epoch (dashed line) and with deep stacked images (continuous 
line). 



not increase the typical crowding of the image. Later on, ar- 
tificial stars spanning small bins of colour and magnitude are 
grouped together to provide estimates of number ratio between 
added and recovered stars - i.e. the completeness - as a func- 
tion of position in the CMDs. 

Typical results of this process are illustrated in the Fig. [TTl 
which shows the completeness as a function of magnitude, for 
both cases of single epoch and deep stacked images. For this 
relatively low density tile, the figure shows that the ~5Q% com- 
pleteness limit goes as deep as 22.20, 21.90, and 21.40 mag in 
the Y, J and bands, respectively. These results are in good 
agreement with the expectations derived from simulated VMC 
images (Kerber et al. I2009I I. 

In order to derive the spatially-resolved SFH of the 
Magellanic Clouds, we are performing a more extensive series 
of artificial star tests in the available images. The results will 
be discussed in detail in a forthcoming paper of this series. 

5.2. RR Lyrae stars and Cepheids 

Radially pulsating stars obey a period-mean density relation 
that forms the basis of their use as standard candles to mea- 
sure distances to the host systems. RR Lyrae stars in particular, 
obey a period-luminosity-metallicity (PLZ) relation in the K 
band which is weakly affected by evolutionary effects, spread 
in stellar mass within the instability strip, and uncertainties in 
reddening corrections (see Longmore et al. ^1986 and Sollima 
et al. 120061 for updates). Similarly, the Cepheid PL relation in 
the K band is much narrower than the corresponding optical 
relations, and less affected by systematic uncertainties in red- 
dening and metal content (Caputo, Marconi & Musella l2000l l. 



The /T-band PLZ relation of RR Lyrae stars has already 
been used by a number of authors to derive the distance to LMC 
clusters and field regions (Dall'Ora et al.'2004', Szewczyk et al. 
|2008, Borissova et al. |2009| ). However, these studies only cover 
a few tiny LMC portions mainly located in or close to the LMC 
bar, and just one LMC cluster, the Reticulum (Dall'Ora et al. 
[2004). In the context of the VMC project, the photometry 
is taken in time series mode in order to obtain mean Ks mag- 
nitudes for RR Lyrae stars and Cepheids over the whole MCs, 
and use their PL relations to measure distances and construct a 
3D map of the entire system. 

Over the last decades, large microlensing surveys such as 
MACHO, EROS and OGLE have been conducted in the op- 
tical domain to search for planetary-to-stellar-mass baryonic 
dark matter in the MW. As a by-product these surveys have 
produced tens of millions of light curves of individual back- 
ground stars. These surveys provide an unprecedented oppor- 
tunity for systematic studies of whole classes of variable stars 
and their host galaxies. They cover different fractions of the 
MC area. Specifically, the MACHO survey (see Alcock et al. 
I2004l and references therein) cover a region mainly encompass- 
ing the LMC and SMC bars. The OGLE surveys (Soszynski et 
al. 120081 120091 and references therein) of which stage IV is in 
progress, cover a progressively larger area extending further 
outside the bar of each Cloud. EROS-2 (Tisserand et al. l2007l l 
is the most extended of these surveys and, at present, is the one 
covering the largest fraction of the VMC survey area. 

We are using the identification, the period and the epoch 
of maximum light of RR Lyrae and Cepheids, identified by 
the microlensing surveys, to fold the /Ts-band light curves pro- 
duced by VMC and derive average Ks magnitudes for these 
variables. In this paper we present results of a preliminary anal- 
ysis based on the combination of VMC, EROS-2 and OGLE 
III data for Cepheids and RR Lyrae stars in two LMC fields, 
namely, the Gaia and 30 Dor fields. 

Given its location on the periphery of the LMC, the Gaia 
field only overlaps with EROS-2. Coordinates, periods and op- 
ticadli ght-curve dataset of RR Lyrae and Cepheids in the Gaia 
field are taken from the EROS-2 catalogue and cross -matched 
to the VMC data. The 30 Dor field is covered by both EROS-2 
and OGLE III, as well as by MACHO, but for the present anal- 
ysis we only employ periods and optical (Johnson-Cousins V, I 
bands) light curves from OGLE III (Soszriski et al. 2008 20()9]l. 

RR Lyrae stars and Cepheids in the Gaia field were ex- 
tracted from the EROS-2 catalogue and the (B - R,R) CMD 
using the following limits: 18.46 < B < 20.03 mag, 0.05 < (B- 
R) < 0.58 mag for the RR Lyrae stars, and 13.39 < B < 17.82 
and 0.89 < P < 15.85 days for the Cepheids. This selection re- 
turned a list of 16337 RR Lyrae stars and 5800 Cepheids over 
the whole field of the LMC covered by EROS-2. RR Lyrae can- 
didates and Cepheids within the Gaia field were then extracted 
according to their coordinates by considering only objects with 
87.8464 < ff < 91.8464 deg and -67.3413 < 6 < -65.3413 
deg, giving 235 RR Lyrae and 47 Cepheid candidates. The 



The EROS blue channel (420-720 nm), overlapping the V and R 
standard bands, and red channel (620-920 nm), roughly matching the 
mean wavelength of the Cousins / band (Tisserand et al. l2007t 



Cioni et al.: The VMC Survey - 1. 



17 



RR Lyrae candidates were then further selected by consider- 
ing only objects with 'proper' periodicities from EROS-2. This 
restricted the sample to 218 sources inside the Gaia field. A 
preliminary cross-match between this catalogue and the Vs43 
was then made using the CrossID query (pairing radius 1"), 
yielding a final catalogue of 117 RR Lyrae candidates and 21 
Cepheids in common. 

Astrometric differences between the EROS-2 sources and 
VMC counterparts is < 1" for 98.4% of the sources (with dif- 
ferences for the remainder in excess of 5")- Periods for these 
stars were checked by analysing their Beros light-curves with 
Graphical Analyser of Time Series (GRATIS) custom soft- 
ware developed at the Bologna Observatory by R Montegriffo 
(Clementini et al. I2000l l confirming the EROS-2 periodicities 
for the majority. Similar methods were employed to select and 
cross-match RR Lyrae stars and Cepheids in the 30 Dor field 
from the OGLE III catalogue, and to extract their band time 
series data from the VMC observations. 

Figure [T2] shows Beros and K^, VMC light curves of a fun- 
damental mode RR Lyrae star (star #15574, with P - 0.601586 
days) and a Cepheid (star #6104, with P - 3.87046 days) in 
the Gaia field. Also shown are Vogle hi and VMC light- 
curves of a fundamental mode RR Lyrae star (star #22926, 
with P = 0.5469449 days) and a Cepheid (star #2871, with 
P = 6.3497921 days) in the 30 Dor fielcS. 

In practise, in the VMC paw-print the same source will (on 
average) appear in at least one more detector within the same 
tile sequence; the magnitudes and uncertainties shown in 
Fig. [12] for RR Lyrae stars are the weighted averages of all 
available observations, so for each night we have one point. 
For the Cepheids all the data available in the VSA release 
(VMC20100607) are shown in the figure, i.e., two points for 
each night, in order to show the excellent raw data for these 
stars. 

Error bars of the individual measurements are shown in 
the figures. The uncertainties for the Cepheids in both fields 
are the same size at the data-points in the figures, while they 
are larger for the 30 Dor RR Lyrae star due to the crowding 
in this part of the LMC. Nevertheless, the typical shape of 
the fundamental mode of the RR Lyrae stars is easily recog- 
nised in Fig. [121 The light-curves are very well sampled for 
both Cepheids and RR Lyrae stars confirming the soundness 
of our observing strategy, and allowing us to derive accurate 
(/Ts) magnitudes without using template light curves. The {Kg) 
magnitudes, computed as the average of the individual mea- 
sures correspond to 18.02 + 0.12 mag, 17.99 + 0.06 mag for 
the Gaia and the 30 Dor RR Lyrae stars, respectively, and to 
14.17 + 0.06 mag, 13.46 ± 0.08 mag for the Gaia and the 30 
Dor Cepheids, respectively where the quoted uncertainties are 
the standard deviations of the averages. 



Database release VMC20 100607 

The variables in the Gaia field have EROS-2 identifications of 
Im0382ml5574 and lm050716104 for the RR Lyrae and Cepheid, re- 
spectively. Similarly, the variables in the 30 Dor field have OGLE 
III identifications of OGLE-LMC-RRLYR-22926 and OGLE-LMC- 
CEP-2871, for the RR Lyrae and Cepheid, respectively. 



The {Ks) magnitudes from the VMC time series will be 
used along with EROS-2 and OGLE III light-curves to con- 
struct PLkZ relations for the RR Lyrae stars and PL, PL-colour 
and Wesenheit relations for Cepheids, to then investigate the 
geometrical distribution of variables in these LMC fields and 
for the MC system in general. 

5.3. Planetary nebulae 

There are approximately 700 objects catalogued as PNe in the 
LMC (Leisy et al. 1997, Reid & Parker 2006a, 2006b) and 
140 in the SMC (Jacoby, private communication; Sanduleak & 
MacConnell & Philip T978^, Jacobv [T980i Jacoby & De Marco 
[2002). Magellanic PNe are best known for their fundamental 
role in the development of the extragalactic standard candle 
[O III] /15007 Planetary Nebula Luminosity Function (PNLF; 
Henize & Westerlund[l963] Jacoby 1980 1989). Distances can 
be measured from the near-universal bright end cut-off across 
all galaxy types, but it remains difficult to explain how old stel- 
lar populations that lack recent star-formation episodes could 
produce progenitors massive enough to power the high cen- 
tral star luminosities populating the bright end (Jacoby 1 19971 
Marigo et al. l2004l Ciardullo 2010). 

Some success in reproducing the observed PNLF has 
been achieved by incorporating the latest hydrodynamic, time- 
dependent models (Schonberner et al. 120071 Mendez et al. 
120081) ■ In these simulations the emphasis is on modelling the 
PNLF rather than explaining it and an alternative explanation 
for [O III] -bright PNe is still required. Ciardullo et al. (120051 ) 
proposed that blue stragglers could evolve into [O III] -bright 
PNe provided they are formed via a close binary interaction. 
Soker (2006) and Frankowski & Soker (2009) suggested sym- 
biotic stars could fulfill the role of [OlII]-bright PNe, but this 
would controversially require the majority of their nebulae to 
be ejected by the white dwarf (see Corradi 120031) . Both of 
these scenarios require alternative binary evolution channels 
for PN formation and as such they are not out of place amongst 
the growing evidence for binarity in PNe (De Marco 2009). 
The frequency of PNe with close binary central stars that went 
through a common-envelope interaction is fairly high at 17 ±5% 
(Miszalski et al. 120091 ) and symbiotic stars could potentially be 
rebranded as wide interacting binary central stars if we insist 
upon the nebula origin as akeady stated. Out of the two sce- 
narios symbiotic stars are more readily accessible despite there 
being a severe paucity of Magellanic symbiotics (Belczyriski et 
al. 120001) . 

Magellanic PNe are well positioned to further advance our 
understanding of the PNLF. Large catalogues of [O III] fluxes 
are available (Jacoby & De Marco 2002 Reid & Parker llOTOl ) 
and the MCs are close enough to allow their PNe to be spatially 
resolved and studied in detail (e.g. Shaw et al. 120061 ). This is a 
critical advantage over more distant [O III] selected populations 
which are frequently assumed to contain only PNe. There are 
however a large variety of potential mimics that can contami- 
nate the PNLF (see Frew & Parker 2010 for a review) and sym- 
biotic stars have been identified in the Local Group (Gongalves 
et al. 2008, Kniazev et al. |2009l ). CiarduUo (l20T0b found a large 
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Fig. 12. Beros and Kg VMC light-curves of the fundamental mode RR Lyrae star #15574 and for Cepheid # 6104 in the Gaia 
SEP field (top) as well as Vogle m and K^^ VMC light-curves of the fundamental mode RR Lyrae star #22926 and for Cepheid # 
2871 in the 30 Dor field (bottom). Data are folded according to the period and epoch of maximum light derived from the EROS-2 
data. 



scatter in the emission line ratios of objects in the top magni- 
tude of the M 33 and LMC Ha PNLFs which supports the case 
that more than one type of object can occupy the bright end. 
The deep near-IR photometry provided by VMC is sensitive to 



the dust associated with many mimics including compact HII 
regions and symbiotic stars. 

The synoptic nature of the VMC survey will detect variabil- 
ity due to Mira pulsations of the most obscured symbiotic stars 
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Fig. 13. Colour-composite images of the LMC PN MG 60 
made from single (left) and stacked (right) exposures of Y 
(blue), J (green) and A^s (red). The effective exposure times 
are 200 sec for all single exposures, and 2400 sec (Y and J) 
and 6000 sec (K^) for the stacked exposures. The images are 
40 X 35 arcsec^ with North up and East to the left. 



that may otherwise be misclassified as PNe. These pulsations 
may not be visible at optical wavelengths (Mikolajewska et al. 
1 19891 1 providing a unique opportunity to increase the number of 
Magellanic symbiotic stars. We might also be sensitive to vari- 
ability in the brightest central stars that could potentially be ev- 
idence for binarity if the variability is later found to be periodic 
(Miszalski et al. 2009), however only a very small fraction of 
the population might be expected to fall within our magnitude 
Umits (Villaver, StanghelUni & Shaw l2007l l. 

Within the first six VMC tiles in the LMC, a combination 
of optical imaging, OGLE light-curves, the VMC near-lR data, 
and SAGE mid-lR observations, reveals that only ~5Q% of the 
98 objects catalogued previously as PNe appear to be genuine. 
These are characterised by the colours 0.4 $ 7 - /Ts ~ 2 and 
0.0 $ y - 7 ^ 0.5. It is encouraging though that almost all 
the genuine PNe are detected in all wavebands and they appear 
especially bright in K^. The strongest emission lines in the near- 
IR for PNe include He I L083yum in Y, Pa /3 in 7, while K, 
contains Br y, multiple He I and molecular H2 lines (Hora et al. 
I1999I I. Figure [13] shows the impact of stacking the individual 
exposures to detect PNe, many of which are only visible in 
the stacked fames. The non-PNe identified in the sample are 
mainly misclassified field stars, compact H II regions or long- 
period variables. Our small fraction of genuine PNe is real and 
not just accounted for by non-detections of faint nebulae in the 
near-IR. It is also preliminary given the statistics are dominated 
heavily by the 30 Dor 6_6 tile whose extensive H II emission 
nebulosity complicates the task of identifying PNe. 

As the VMC survey becomes more complete, near-IR lumi- 
nosity functions will be constructed and compared with their 
[OIII] (Jacoby & De Marco l2002l Reid & Parker l 20T0] l and 
mid-IR (Hora et al. I2008I I counterparts. Improved diagnostic 
capabilities with the sensitivity of the deep stacks and SAGE 
mid-IR photometry will also be applied to discover new PNe 
inside and especially outside the central 5x5 degrees cov- 
ered by Reid & Parker to create a more complete census of 
Magellanic PNe with accurate coordinates. 



5.4. The Stellar cluster KMHK 1577 

Stellar clusters are among the primary targets of the VMC sur- 
vey. The detection of known stellar clusters will be examined, 
and new clusters will be searched for using a method simi- 
lar to that already adopted in the study of the MW (Ivanov et 
al. 120021 Borissova et al. I2003I I. The analysis of stellar clus- 
ters will be centred on the study of CMDs to estimate their 
ages and metallicities as well as on the comparison with results 
obtained from optical surveys: MCPS (Zaritsky et al. I2004I I. 
OGLE (Pietrzynski et al. [T9981 [19991 ). those by de Grijs & 
Anders ("20061 and other dedicated studies. Ages, masses and 
metallicities of stellar clusters will allows us to discuss the SFH 
of the MCs (e.g. synchronised bursts at ~ 0.2 and 2 Gyr be- 
tween the LMC and the SMC) and radial abundance gradients. 
The first step of this work has been the compilation of a cata- 
logue of known stellar clusters located in the Gaia SEP and 30 
Dor fields from the Hst of Bica et al. (I2008I I. 

In this paper we show preliminary results obtained from the 
study of one cluster in the 8_8 field: KHMK 1577. This cluster 
appears in the literature in different catalogues, including that 
from Bica et al. (2008), but it is poorly studied and its properties 
(i.e. age and metallicity) are unknown. It was chosen because of 
its favourable location at the centre of a VMC paw-print which 
maximises the availability of data for its stellar members. The 
VMC observations, available in this field, are listed in Tabs.[4| 
and[5] The CMD of the cluster and its immediate surroundings 
were examined using photometry from deep VSA stacks. The 
three panels in Fig. [14] show images of the cluster region in 
the Y, J and bands where the elliptical region occupied by 
KMHK 1577 is indicated. 

The sources in the cluster region were divided into two con- 
centric groups: an inner group comprising all sources within an 
ellipse with a major axis of 0.75', a minor axis of 0.70' and 
a position angle of the major axis of 160° (Bica et al. I2008I I. 
and an outer group of sources within a circular area (0.63' in 
radius) equal in size to twice the area of the inner ellipse. The 
stars within the elliptical region are likely members of the clus- 
ter, although confirmation from radial velocities would be nec- 
essary, while stars in the outer region are unlikely associated 
with the cluster. By extracting the VMC data from the VSA we 
find 69 sources in the inner group and 98 sources in the outer 
group. If only sources detected in all three wave bands are con- 
sidered then there are 41 and 55 sources in the inner and outer 
group, respectively. The inner group has a higher source den- 
sity. Using the VSA source classification flag there are 21 stars 
and 20 galaxies in the inner group and 34 stars and 21 galaxies 
in the outer group. Extended objects at faint magnitudes, how- 
ever, are most probably too faint for the classification criterion 
into stars and galaxies to work. These objects may be elon- 
gated either because of a low signal-to-noise or because they 
are blends. Their nature will be clearer once all VMC epochs 
have been obtained and stacked. Only one very bright object 
is classified as extended and this may be the result of blending 
in the cluster centre. Objects with red colours and low lumi- 
nosities are likely extended objects (galaxies, PNe and YSOs), 
none are present inside the ellipse. Sources in the outer region. 
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Fig. 14. KMHK 1577 observed by VISTA in the Y, J and bands as part of the VMC survey: these images are about 2' x 2' in 
size and refer to 1200, 1400 and ~ 4000 sec of integration, respectively. The ellipse shows the size of the cluster as indicated by 
Bica et al. dlOOSl) . 
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Fig. 15. Magnitude differences between 2MASS and VMC de- 
tections in the J (empty squares) and (filled triangles) bands. 



that are homogeneously distributed across the area, are either 
giant stars of the LMC field population or extended objects. 

A comparison with 2MASS sources detected within the el- 
liptical region defining the cluster shows that at bright magni- 
tudes (K^ ~ 13) there is a very good agreement between the 
2MASS and VMC magnitudes, with increasingly discrepant 
values when approaching the 2MASS sensitivity limit. For the 
sources in common between VMC and 2MASS (Fig. [TsT l the 
magnitude difference in J and Ks varies between being practi- 
cally negligible to values of 0.4 - 0.5 mag. The difference in 
the J - colour is up to 0.8 mag. 

The CMD shown in Fig. [T6] indicates that the present sen- 
sitivity Umit of VMC is at A's ~ 21 mag; this value agrees 
with expectations. Figure[T6]shows that the entire RGB is well 
detected as well as the sub-giant branch and stars at the MS 
turn-off point. Photometric uncertainties vary with brightness 
and in the Ks band sources with a magnitude of 19 - 20 have 
uncertainties of 0.2 - 0.4 mag, sources with \6 < < 17 
have uncertainties of 0.01 - 0.04 mag and the brightest sources 
at /Ts ~ 13 mag have an uncertainty of ~ 0.003 mag. Y and 
J bands sources have smaller photometric uncertainties than 
sources in K^. 
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Fig. 16. CMD (J - Ks, Ks) of the inner region around the stel- 
lar cluster KMHK 1577. Points represent stellar (filled trian- 
gles) and non-stellar (asterisks) objects. The line is the best fit 
isochrones from Marigo et al. (,2008j with an age of 0.63 Gyr 
and a metallicity of Z = 0.003. 

Using theoretical isochrones by Marigo et al. (120081 ) an 
age for the cluster of 0.63 + 0.10 Gyr and a metallicity of 
Z= 0.003 + 0.001 were estimated. These values were derived 
from the average of the best fit isochrones found in three dif- 
ferent CMDs: (J - Ks, Ks), (Y - Ks, Ks) and (Y - J, J). In 
each of these the best fit age was the same but the metallic- 
ity varied within the given dispersion. A similar procedure ap- 
plied to the outer region gives an older (1.4 Gyr) and more 
metal poor (Z- 0.0004) population. The absorption due to 
both foreground and interstellar extinction was estimated from 
Ay - 0.55 mag (Zaritsky et al. 2004) and the A^/Ay ratios for 
an 'average' LMC according to Gordon et al. (I2003I I. This re- 
sults in Ay - 0.22, Ay = 0.14 and Ak, - 0.02 mag where the 
Y band value is an extrapolation from the nearest bands. The 
age of the cluster is typical of young LMC clusters while the 
low metallicity agrees with the cluster being located in a region 
quite far from the LMC centre. 

Figure [17] shows a CMD (B - V, B) for the cluster KHMK 
1577 from the MCPS data (Zaritsky et al. l2004l l. With refer- 
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Fig. 17. CMD (y, B-V) obtained from the MCPS data 
(Zaritsky et al. l2004l l. Points, for VMC counterparts, are as in 
Fig. Hi 

ence to the two groups analysed here we find 63 sources in 
the inner group and 70 sources in the outer group. Ahhough it 
appears that MCPS has more stars in the central region of the 
cluster and allows for a better identification of the density en- 
hancement, the VMC data delineate better features in the CMD 
(compare Fig.[T6]with Fig.fTTTi. Extinction and crowding effects 
in the optical data are the most likely explanation for the differ- 
ence between the two diagrams. 

The VMC data harbour a great potential for investigating 
the evolutionary properties of stellar clusters. The internal spa- 
tial distribution of stars for each individual stellar cluster would 
provide information for dynamical studies. 



The VMC survey will be of immense value to the astro- 
nomical community because the data will represent the only 
counterpart for existing optical surveys at a similar sensitivity 
(e.g. MCPS) and for the large number of unclassified objects 
observed by the Spitzer Space Telescope in the mid-IR (Blum 
et al. 2006 ). Note that the near-IR 2MASS survey has observed 
only about 6% of the stars that the VMC survey is expected to 
detect. The VMC data cover the bulk of the Magellanic system, 
as opposed to the tiny regions sampled by the Hubble Space 
Telescope, and the limited area covered by most of the other 
dedicated, ground-based observations at the same sensitivity. 

Among the other VISTA Public Surveys, the VISTA 
Hemisphere Survey (VHS) will also contribute to the investiga- 
tion of the Magellanic system. The VHS is ~3 mag shallower 
than the VMC survey but it covers an extended area around the 
Magellanic system to complement the VMC data. 
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and beyond' (P.I. M. Marconi). 



6. Conclusions 

The VMC survey is a homogeneous and uniform YJKf^ sur- 
vey of ~184 deg^ across the Magellanic system (Fig. |2]). It is 
an ESO Public Survey that started observations in November 
2009 and will run for approximately five years. The VMC sur- 
vey parameters are described in Tab. [1] Images and catalogues 
will be delivered to the astronomical community at regular in- 
tervals with the first release expected in 2011. The VMC data 
will provide, among other things, a detailed history of star for- 
mation across the Magellanic system and a measurement of its 
3D geometry. 

This paper presents the VMC survey strategy and first re- 
sults aimed at assessing their scientific quality. These show the 
potential of the survey in addressing its main science goals and 
validates the expected sensitivity of the VMC data. Colour- 
magnitude and colour-colour diagrams show a wealth of sub- 
structures and a clear separation from Galactic foreground 
stars. These diagrams will form the core of the SFH analysis. 
To illustrate some of the scientific applications of the VMC sur- 
vey, Cepheids and RR Lyrae stars are shown to display clear 
near-IR light-curves, PNe are detected thanks to the deep 
band images dominated by Bry emission, and stellar cluster 
parameters are derived from best fit isochrones to the VMC 
colour-magnitude diagrams. 
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Appendix A: VMC tile centres 

Tables lA.ll and IA.2I show the centres of VMC tiles covering 
the LMC, while Table IA.3I shows the centres of tiles covering 
the SMC and Table I A.4l shows the centres of tiles covering the 
Bridge. The tile identification is formed of two numbers. The 
first number indicates the row and the second number the col- 
umn that correspond to the location of a given tile. Row num- 
bers increase from bottom to top while column numbers in- 
crease from right to left. Refer to Figures lAniA.2l and lA.3l for 
the location of tiles across the LMC, SMC and Bridge compo- 
nents of the Magellanic system, respectively. For the Stream the 
central coordinates of two distinct tiles are indicated in Table 
IA.5I Note that since tiles were generated from almost rectan- 
gular grids, some outer tiles were removed leading to the tile 
numbers shown here. 



Table A.l. LMC tile centres 



Tile 




a 




6 




Comments 


2.3 


04 


48 


04.752 


-74 


54 


11.880 




2.4 


05 


04 


42.696 


-75 


04 


45.120 




2.5 


05 


21 


38.664 


-75 


10 


50.160 




2.6 


05 


38 


43.056 


-75 


12 


21.240 




2.7 


05 


55 


45.720 


-75 


09 


17.280 




3.2 


04 


37 


05.256 


-73 


14 


30.120 




3.3 


04 


51 


59.640 


-73 


28 


09.120 




3.4 


05 


07 


14.472 


-73 


37 


49.800 




3.5 


05 


22 


43.056 


-73 


43 


25.320 


started 


3.6 


05 


38 


18.096 


-73 


44 


51.000 




3.7 


05 


53 


51.912 


-73 


42 


05.760 




3.8 


06 


09 


16.920 


-73 


35 


12.120 




4.2 


04 


41 


30.768 


-71 


49 


16.320 


started 


4.3 


04 


55 


19.512 


-72 


01 


53.400 


started 


4.4 


05 


09 


24.288 


-72 


10 


49.800 




4.5 


05 


23 


39.816 


-72 


15 


59.760 




4.6 


05 


38 


00.408 


-72 


17 


20.040 


started 


4.7 


05 


52 


20.064 


-72 


14 


49.920 




4.8 


06 


06 


32.952 


-72 


08 


31.200 




4.9 


06 


20 


33.408 


-71 


58 


27.120 




5.1 


04 


32 


43.848 


-70 


08 


40.200 




5.2 


04 


45 


19.440 


-70 


23 


43.800 




5.3 


04 


58 


11.664 


-70 


35 


27.960 




5.4 


05 


11 


16.704 


-70 


43 


46.200 




5.5 


05 


24 


30.336 


-70 


48 


34.200 


started 


5.6 


05 


37 


47.952 


-70 


49 


49.440 




5.7 


05 


51 


04.872 


-70 


47 


31.200 




5.8 


06 


04 


16.416 


-70 


41 


40.560 




5.9 


06 


17 


18.096 


-70 


32 


20.760 




6.1 


04 


36 


49.488 


-68 


43 


50.880 




6.2 


04 


48 


39.072 


-68 


57 


56.520 




6.3 


05 


00 


42.216 


-69 


08 


54.240 




6.4 


05 


12 


55.800 


-69 


16 


39.360 


started 


6.5 


05 


25 


16.271 


-69 


21 


08.280 




6.6 


05 


37 


40.008 


-69 


22 


18.120 


30 Dor - completed 


6.7 


05 


50 


03.168 


-69 


20 


09.240 




6.8 


06 


02 


21.984 


-69 


14 


42.360 


in queue 


6.9 


06 


14 


32.832 


-69 


05 


59.640 




6.10 


06 


26 


32.280 


-68 


54 


05.760 




7.2 


04 


51 


34.992 


-67 


31 


57.000 




7.3 


05 


02 


55.200 


-67 


42 


14.760 


in queue 


7.4 


05 


14 


23.976 


-67 


49 


30.720 




7.5 


05 


25 


58.440 


-67 


53 


42.000 


in queue 


7.6 


05 


37 


35.544 


-67 


54 


47.160 




7.7 


05 


49 


12.192 


-67 


52 


45.480 




7.8 


06 


00 


45.240 


-67 


47 


38.040 




7.9 


06 


12 


11.736 


-67 


39 


26.640 




7.10 


06 


23 


28.800 


-67 


28 


14.880 
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Table A.2. LMC tile centres (continue) 



Tile a 6 Comments 



8.2 


04 


54 


11.568 


-66 


05 


47.760 




8.3 


05 


04 


53.952 


-66 


15 


29.880 


started 


8.4 


05 


15 


43.464 


-66 


22 


19.920 




8.5 


05 


26 


37.704 


-66 


26 


15.720 




8.6 


05 


37 


34.104 


-66 


27 


15.840 




8.7 


05 


48 


30.120 


-66 


25 


19.920 




8.8 


05 


59 


23.136 


-66 


20 


28.680 


Gaia - completed 


8.9 


06 


10 


10.632 


-66 


12 


43.560 




9.3 


05 


06 


40.632 


-64 


48 


40.320 


started 


9.4 


05 


16 


55.464 


-64 


55 


07.680 




9.5 


05 


27 


14.256 


-64 


58 


49.440 




9.6 


05 


37 


34.872 


-64 


59 


44.520 




9.7 


05 


47 


55.128 


-64 


57 


52.920 


in queue 


9.8 


05 


58 


12.816 


-64 


53 


15.000 




9.9 


06 


08 


25.848 


-64 


45 


52.560 




10.4 


05 


18 


01.536 


-63 


27 


54.000 




10.5 


05 


27 


48.912 


-63 


31 


22.800 




10.6 


05 


37 


37.800 


-63 


32 


13.200 




10.7 


05 


47 


26.352 


-63 


30 


24.840 




11.6 


05 


37 


42.432 


-62 


04 


41.520 





Table A.3. SMC tile centres 



Tile 




a 




5 




Comments 


2.2 


00 


21 : 


43.920 


-75 


12 


04.320 




2.3 


00 


44 : 


35.904 


-75 


18 


13.320 




2.4 


01 


07 : 


33.864 


-75 


15 


59.760 




2.5 


01 


30 : 


12.624 


-75 


05 


27.600 




3.1 


00 


02 : 


39.912 


-73 


53 


31.920 




3.2 


00 


23 : 


35.544 


-74 


06 


57.240 




3.3 


00 


44 : 


55.896 


-74 


12 


42.120 


started 


3.4 


01 


06 : 


21.120 


-74 


10 


38.640 




3.5 


01 


27 : 


30.816 


-74 


00 


49.320 


started 


3.6 


01 


48 : 


06.120 


-73 


43 


28.200 




4.1 


00 


05 : 


33.864 


-72 


49 


12.000 




4.2 


00 


25 : 


14.088 


-73 


01 


47.640 




4.3 


00 


45 : 


14.688 


-73 


07 


11.280 




4.4 


01 


05 : 


19.272 


-73 


05 


15.360 


NGC419 


4.5 


01 


25 : 


11.088 


-72 


56 


02.760 




4.6 


01 


44 : 


34.512 


-72 


39 


44.640 




5.2 


00 


26 : 


41.688 


-71 


56 


35.880 


started 


5.3 


00 


45 : 


32.232 


-72 


01 


40.080 




5.4 


01 


04 : 


26.112 


-71 


59 


51.000 


NGC411 - 


5.5 


01 


23 : 


09.336 


-71 


51 


09.720 




5.6 


01 


41 : 


28.800 


-71 


35 


47.040 


in queue 


6.2 


00 


28 : 


00.192 


-70 


51 


21.960 




6.3 


00 


45 : 


48.792 


-70 


56 


09.240 




6.4 


01 


03 : 


40.152 


-70 


54 


25.200 




6.5 


01 


21 : 


22.560 


-70 


46 


11.640 




7.3 


00 


46 : 


04.728 


-69 


50 


38.040 




7.4 


01 


03 : 


00.480 


-69 


48 


58.320 





Table A.4. Bridge tile centres 



Tile 




a 




6 




Comments 


1.2 


01 


49 : 


51.960 


-74 


43 : 


31.800 




1.3 


02 


11 : 


35.232 


-75 


05 : 


04.560 




2.3 


02 


14 : 


46.584 


-74 


00 : 


47.520 


started 


2.4 


02 


35 : 


28.440 


-74 


13 : 


18.840 




2.7 


03 


39 : 


43.800 


-74 


04 : 


51.960 




2.8 


04 


00 : 


21.072 


-73 


46 : 


37.560 


started 


2.9 


04 


20 : 


05.640 


-73 


21 : 


14.040 




3.3 


02 


17 : 


35.496 


-72 


56 : 


22.200 




3.4 


02 


37 : 


24.888 


-73 


08 : 


16.440 




3.5 


02 


57 : 


33.288 


-73 


12 : 


52.200 


started 


3.6 


03 


17 : 


43.776 


-73 


10 : 


03.720 




3.7 


03 


37 : 


39.240 


-72 


59 : 


54.600 


in queue 


3.8 


03 


57 : 


03.888 


-72 


42 : 


37.800 





Table A.5. Stream tile centres 



Tile 




a 






S 




Comments 


1.1 


03 


30 : 


03.87 


-64 


25 


24.8 


started 


2.1 


00 


11 : 


59.30 


-64 


39 


31.8 


started 




7.5h 7h 6.5h Bh 5.5h 5h 4.5^ 4h i.5h 



a (J2000) 

Fig.A.l. LMC area tiled for VMC observations. The under- 
lying image shows the Ha distribution (Gaustad et al. 1200 11 1. 
VISTA tiles are colour-coded as in Fig. |2] 
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2,5h 2h l.5h Ih O.Sh Oh 23.5h 

a (J2000) 

Fig.A.2. SMC area tiled for VMC observations. The under- 
lying image shows the Ha distribution (Gaustad et al. 1200 11 1. 
VISTA tiles are colour-coded as in Fig.|2] 
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Fig. A.3. Bridge area tiled for VMC observations. The underly- 
ing image shows the Ho- distribution (McClure-GrifRths et al. 
|2009] l. VISTA tiles are colour-coded as in Fig. |2] 



